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Abstract 



A search is performed for an excess of events, over the standard model expecta- 
tions, with a photon, a lepton, and large missing transverse energy in pp collisions 
at = 7 TeV. Such events are expected in many new physics models, in particular 
a supersymmetric theory that is broken via a gauge-mediated mechanism, when the 
lightest charged and neutral gauginos are mass degenerate. The data sample used in 
this search corresponds to an integrated luminosity of 35 pb^^ collected with the CMS 
detector at the LHC. No evidence of such an excess above the standard model back- 
grounds, dominated by W7 production, is found. The results are presented as 95% 
confidence level upper limits on the cross section for a benchmark gauge-mediated 
scenario, and are then converted into exclusion limits on the squark, gluino, and wino 
masses. 
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1 Introduction 

Supersymmetry (SUSY) is one of the best studied scenarios for physics beyond the standard 
model (SM). It alleviates the hierarchy problem and may provide a path to the grand unifica- 
tion of forces. Of particular interest is a scenario with gauge-mediated SUSY breaking (T-'9], 
which elegantly addresses the SUSY flavour problem. In the gauge mediation scenario, SUSY 
breaking occurs at energy scales much smaller than the Planck scale, resulting in the gravitino 
(G) as the lightest SUSY particle (LSP), which would escape detection and result in an imbal- 
ance of momentum in the plane transverse to the beam direction (E™*®). In one scenario of 
general gauge-mediated (GGM) SUSY \M^, the next-to-lightest SUSY particles (NLSP) are 
mostly the superpartners of the SU{2)i gauge fields (winos). This scenario results in a small 
mass splitting between the charged (W^) and neutral (W^) winos. The charged wino decays 
promptly into a gravitino and a W boson, and the neutral wino into a gravitino and a photon or 
a Z boson [12J. If R-parity is conserved, coloured SUSY particles are pair-produced, resulting 
in the presence of two NLSPs and multiple jets per event. An example process is shown in 
Fig. [l] The event topologies of interest therefore contain, in addition to EJ^'^®, one of the follow- 
ing signatures: two photons (7); lepton(s) {£) and one photon; lepton(s) and jets; one photon 
and jets; or all jets. Final states with two or more leptons are suppressed because of the W and 
Z leptonic branching fractions. 

In this paper, we perform the first search for SUSY in the t + j + Ei^'^® + X final state produced 
in pp collisions at = 7TeV at the Large Hadron Collider (LHC). The i denotes an electron 
or a muon. Previous searches for anomalous production of such events were performed at the 
Fermilab Tevatron IITSHTSi . 




Figure 1: Feynman diagram for an example of a process with strong production of charged 
and neutral winos (W^ and W^) as co-NLSP with multiple jets (j), resulting in the signature 

e + j + E^ss + X. 



2 The CMS Detector 

The Compact Muon Solenoid (CMS) is a nearly hermetic, multipurpose detector and is de- 
scribed in detail elsewhere [16|. Charged particle trajectories are measured with silicon pixel 
and strip tracking detectors encompassed by a 3.8 T superconducting solenoid, with full az- 
imuthal coverage and < 2.5. The azimuthal angle (p is defined in the plane transverse to the 
beam direction. The pseudorapidity t] is defined as t] = — ln[tan(0/2)], where 6 is the polar 
angle with respect to the counterclockwise proton beam direction, with the origin at the center 
of the detector. Energy deposits of particles are measured with lead-tungstate crystal electro- 
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4 Background Estimation 



magnetic (ECAL) and brass /scintillator (\f]\ < 3) or iron /quartz-fiber (3 < l^yl < 5) hadronic 
(HCAL) calorimeters. Muon candidates are identified in gas-ionisation detectors embedded in 
the steel return yoke of the magnet. 

3 Candidate Event Selection 

We use the data sample collected by the CMS detector at the LHC from March to November 
2010 corresponding to an integrated luminosity of 35 pb^^ [17J. Events are required to have at 
least one high transverse momentum (pj) electron or muon, in addition to a photon. Events 
are selected and recorded using a two-tiered trigger system. We use a suite of single-electron 
and single-muon triggers with several thresholds changing to match the rapidly increasing in- 
stantaneous luminosity of the LHC. The electron candidates are reconstructed from clusters of 
energy in the ECAL barrel (|f/| < 1.479) and endcaps (1.479 < \r]\ < 3) and are required to 
match charged tracks reconstructed in the silicon tracking detectors. Muons are reconstructed 
as charged tracks matched to hit patterns in the muon detectors. We further require the lead- 
ing electron or muon to have px > 20GeV, \f]\ < 2.1 and to satisfy a set of identification and 
isolation criteria [18J . For electrons, the region from 1.44 < |7/| < 1.57 is excluded since it 
contains services and cables exiting between the transition of the barrel and endcap calorime- 
ters and leads to lower quality reconstructed clusters. Events in which the electron is con- 
sistent with having originated from a photon conversion are rejected. The efficiency of the 
triggers for events with an electron (muon) with pj above 20 GeV is measured to be 0.98 ± 0.01 
(0.968 ± 0.004). Lepton identification and reconstruction efficiencies are measured in data using 
Z — )■ £^i^ decays. The ratio of the lepton efficiencies measured in data to those from Monte 
Carlo (MC) simulation (described below) is found to be 0.928 ± 0.015 (0.990 ± 0.001) for an 
electron (muon). We require the presence of at least one photon with pj > 30 GeV, |?/| < 1.44, 
spatially separated from the lepton by AR(£, 7) = {Arj)^ + (A^)^ > 0.4. The reconstruction 
of photon candidates in the calorimeter is similar to that of the electrons, with the additional 
requirement that they do not have associated hit patterns in the pixel detector consistent with 
a track. The data/MC scale factor for the photon efficiency is 0.97 ± 0.02. We demand that the 
lepton and the photon candidates be isolated from any other activity in the tracker, ECAL, or 
HCAL. The events are required to contain at least one reconstructed primary vertex within 24 
cm of the geometric center of the detector along the beamline and a transverse distance with 
respect to the reconstructed beam axis of less than 2 cm. These selections define our candidate 
sample, which consists of 264 (126) events in the e'y (jiy) channel. Finally, in order to minimise 
the number of expected background events while retaining high signal efficiency, we select 
events with £^^^^ > 100 GeV. The £^'®* is reconstructed using calorimeter energy deposits and 
is corrected for the contribution from the muons and for the expected calorimeter response for 
each reconstructed hadron track |[19l . 

4 Background Estimation 

The dominant background process is W7 production, with the W boson decaying leptonically. 
This background is estimated using the madgraph MC generator ||20| with parton showering 
using PYTHIA [21J. The result is scaled from the leading-order (LO) to the next-to-leading-order 
(NLO) cross section obtained from the Baur NLO W7 generator [22J and using the CTEQ66 
NLO parton distribution function (PDF) sets [23 J. For example, for a photon with pj = 50 GeV 
the NLO/LO K factor is 2.3. The NLO calculation is in good agreement with the W7 produc- 
tion cross section measured by the CMS collaboration [2^, over the relevant phase space of 
that measurement, giving confidence in the NLO prediction used here. The MC events are 
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processed with the full CMS detector simulation based on the Geant4 p5| package. To esti- 
mate the uncertainty on the K factor due to higher-order corrections, we vary the factorisation 
and renormalisation scales by a factor of two and observe a variation in the cross section of 
approximately 20%. Additional systematic uncertainties on this background estimate include 
the uncertainties on the trigger efficiencies and data/MC scale factors described above, the 
PDF set estimated using the PDF4LHC working group procedure (2%) ||261 , and the integrated 
luminosity measurement (4%) |IT7| . 

The next most important instrumental background originates from SM processes that have 
misidentified leptons and photons. Data-driven methods are used to determine all backgrounds 
involving object misidentification Il27t . We consider sources of background in which a jet or an 
electron is misidentified as a photon (jet — > 7 and e — > 7) or in which a jet is misidentified 
as a lepton. The dominant SM processes that contribute to the jet — > 7 misidentification back- 
ground include W+jets and QCD multijet production. The jet — > 7 background is determined 
by selecting a sample of events that satisfy the criteria of our £7 candidate sample except that 
the photon is required to fail either the isolation or electromagnetic shower shape criterion. We 
then weight this sample by the probability of misidentifying a jet as a photon [28 1 to estimate 
the total contribution of this background to the signal sample. The SM processes contributing 
to the e — 7- 7 misidentification background are predominantly Z and tt production. The e — )■ 7 
background is determined in a similar way to the jet — ?► 7 background, by selecting a sample 
of events that satisfy the criteria of our £7 candidate sample but with the additional require- 
ment that the photon be matched to a pattern of hits in the pixel tracker. We then weight this 
sample by the probability of misidentifying an electron as a photon (fjj to estimate the total 
contribution of this background to the signal sample. The probability, fy^, is determined from 
Z e+e data to be 0.014 ± 0.004. 

The strategy for determining QCD contributions which have poorly measured £™^^^ is to se- 
lect events with kinematic properties similar to the candidate events, but which are known to 
have no genuine £^^^^. A control sample of events dominated by Z — )■ e+e^ decays is used. 
These events are further reweighted to reproduce the £7 sample kinematics. In particular, we 
reproduce the £7 transverse energy distribution and, therefore, the transverse energy of the 
hadronic recoil against the £7 candidates, ensuring a correct description of the £™'®®. In order 
to describe the distribution of the transverse mass of the lepton and £^'®® system, defined as 

Mt = ^2 ■ pt/, ■ £™'®® • (1 — cos{(p( — ^grmss)), the lepton pj distribution is also reweighted. We 

determine the background contribution by normalising the control sample in such a way that 
the total background matches the number of selected events in data with £^'** less than 30 GeV. 
Tables [l] and |2] summarise the breakdown of the backgrounds for E^^^^ > 40 GeV and for our 
optimal selection of £^^^^ > 100 GeV. The dominant sources of systematic uncertainties are due 
to the limited precision of the data-driven methods and the uncertainties on the corresponding 
misidentification rate measurements. 

5 Results 

The £^'®® distribution is shown in Fig. |2] (top) for the e7 and fij channels combined, together 
with the sum of the expected background contributions, as well as the expectation of the signal 
from a SUSY benchmark point (referred to as GMC, described below) for comparison. Figure |2] 
(bottom) shows the Mj distribution. The data agree well with the SM expectations and do not 
exhibit any evidence for an excess of events with high £^'^®. 

We performed several cross-checks to support our modeling of the SM backgrounds. For exam- 
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5 Results 



Table 1: Event counts for data and expected background in the ej channel. The uncertainty on 
the total background includes the anti-correlation between the uncertainties on the individual 
components due to the scaling of the QCD background to the events with E^®^ < 30 GeV 
as described in the text. The expected event yields from the SUSY benchmark point GMC 
(Section ISll are also shown. 





No E^'^^ selection 


E^'ss > 40 GeV 


E^i'^s > 100 GeV 




44.5 ± 9.2 


16.1 ±3.4 


1.68 ±0.42 


jet^ 7 


20.3 ± 4.5 


3.1 ±0.9 


0.02 ± 0.02 


e — > 7 


70.5 ± 19.1 


0.3 ±0.1 


0.04 ± 0.03 


QCD 


134 ± 28 


0.4 ±0.2 


0.00 ± 0.00 


Total background 


269 ± 18 


19.9 ±3.7 


1.74 ±0.43 


data 


264 


16 


1 


SUSY GMC prediction 


3.94 ± 0.79 


3.76 ± 0.75 


2.79 ± 0.56 



Table 2: Event counts for data and expected background in the jij channel. The uncertainty on 
the total background includes the anti-correlation between the uncertainties on the individual 
components due to the scaling of the QCD background to the events with EJ^'^^ < 30 GeV 
as described in the text. The expected event yields from the SUSY benchmark point GMC 
(Section ISll are also shown. 





No E^^^« selection 


E™«« > 40 GeV 


E^'^s > 100 GeV 


W7 


44.8 ± 9.3 


15.9 ± 3.4 


1.40 ±0.37 


jet^ 7 


18.0 ±4.0 


3.7 ±1.1 


0.10 ±0.09 


e — > 7 


1.2 ±0.4 


0.6 ±0.2 


0.09 ± 0.04 


QCD 


58.3 ± 15.1 


0.2 ±0.1 


0.00 ± 0.00 


Total background 


122.3 ± 12.3 


20.4 ± 3.7 


1.59 ±0.39 


data 


126 


27 


1 


SUSY GMC prediction 


5.12 ± 1.02 


4.84 ± 0.96 


3.66 ± 0.73 
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pie, events with Mj above 70 GeV have small QCD background, as shown in Fig. |2] (bottom), 
and we verified that we can describe all the kinematic distributions in these events with the 
sum of the expected backgrounds, except the contribution from QCD. We also checked that the 
photon candidates have a timing consistent with having originated from the primary vertex, to 
verify that non-collision backgrounds are negligible. Finally, we verified that our QCD back- 
ground determination procedure works in the presence of a large range of SUSY signals by 
mixing simulated SUSY events into the data sample and recalculating the QCD background; 
the change in its value was found to be negligible. 

Since we observe no evidence of an excess of events with high compared to the SM ex- 
pectations, we interpret our results by placing 95% confidence level (CL) lower limits on GGM 
SUSY particle masses in a benchmark model. The GGM SUSY [12] signal is simulated using 
the PYTHIA [21J MC generator version 6.420 with tune D6T [29J, and assuming that the squark 
mass im{q)) is approximately equal to the gluino mass im{g)), with the charged and neutral 
winos as co-NLSPs, and tan j6 = 2. This search focuses on the m{W^) > 100 GeV region, which 
is beyond the sensitivity of the previous collider experiments llT3HT5ll30| . We therefore generate 
points with m{g) and m{q) in the range of 350-800 GeV and ?n(W^) ~ jn(W^) in the range of 
100^80 GeV. The generated events are processed with the full CMS detector simulation using 
the Geant4 package. The NLO K factors were calculated using PROSPINO [31 J and are found on 
average to be 1.4 with a ~ 20% variation. As an example of a SUSY benchmark model, we use 
the point generated with m{g) = m{q) = 450 GeV and m(W'') m(W^) = 195GeV, referred 
to as CMC. 

We use a Bayesian approach to determine the 95% CL upper limits on the cross sections in the 
benchmark SUSY model described above, combining the likelihoods of the e'y and f/7 channels. 
The inputs for both channels are the sum of the data-driven backgrounds, the W7 background 
as determined from the MC simulation and corrected by the appropriate data/MC scale fac- 
tors, and, for each MC SUSY point, the efficiency times acceptance, corrected for the appropri- 
ate data/MC scale factors. Each of these inputs for a given channel has log-normal distributed 
nuisance parameters that vary by the uncertainty in each quantity. We consider the 4% uncer- 
tainty on the luminosity as an additional nuisance parameter. As an example of the calculation, 
the efficiency times acceptance times branching ratio for the GMC benchmark point for ey (jiy) 
is 0.0063 ± 0.0003 (0.0083 ± 0.0004). The 95% CL upper limit on the cross section obtained for 
this point is 8.7 pb. The limit is less than the predicted NLO cross section for this point, (12.7 ± 
2.6) pb, with uncertainties due to the PDF [26J (10%) and the renormalization scale (18%). We 
therefore exclude this point with more than 95% confidence. The observed 95% CL upper limits 
on the cross section are shown in Fig.|3]as a function of the squark / gluino and the wino masses 
for the kinematic selections applied in this analysis (an electron with pj > 20 GeV, \rj\ < 1.44 
or 1.57 < \r]\ < 2.1; a muon with pj > 20GeV and \r]\ < 2.1; a photon with pj > 30GeV, 
\rj\ < 1.44; events with E^^^^ > 100 GeV). Figure |4] displays the resulting 95% CL exclusion re- 
gion in the squark / gluino and wino mass plane, where masses in the area below the curve are 
excluded. The shaded region around the central value of the exclusion contour represents the 
variation of the predicted cross section by one standard deviation. Also shown is the expected 
95% CL exclusion as well as the the region excluded by previous searches at the Tevatron in the 
same channel |[T3tiT5| , as interpreted in ||30|. 

6 Conclusions 

We performed the first search in proton-proton collisions at ^/s = 7 TeV for an excess of events 
with a photon, a lepton, and large E^^ over the SM expectation. Using a data sample corre- 
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6 Conclusions 




20 40 60 80 100 120 140 160 180 200 



lA^ (GeV) 

Figure 2: The EJ^'^® (top) and Mj distributions (bottom) for the combined and }ij samples, 
with logarithmic scales on the vertical axes; the last bin shown in both distributions includes 
the overflow events. The black points represent the 35 pb^^ data sample, and the coloured 
histograms show the individual background components. The total background and its uncer- 
tainty are represented by the shaded band. The expectation from the SUSY benchmark point 
GMC is shown by the red dashed line. 
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wino mass (GeV) 

Figure 3: The observed 95% CL upper limits on the cross section as a function of the 
squark / gluino mass (since m{q) ^ ni{g)) and the wino mass is shown for the range of MC 
points generated for this analysis. The star indicates the SUSY benchmark point GMC. 
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> 

O 

(0 

E 

o 
c 



(0 
3 

0) 




400 500 600 700 

wino mass (GeV) 

Figure 4: The 95% CL exclusion region as a function of the squark / gluino mass (since m{q) 
tn[g)) and the wino mass. The area below the curve is excluded by this analysis. The shaded 
region around the central value of the exclusion contour reflects the uncertainty in the NLO 
predicted cross section. The dashed line represents the expected 95% CL limit from the MC 
simulation. The region excluded by the Tevatron in this channel l|T3l - [T5l as interpreted in t30] 
is also shown. The star indicates the SUSY benchmark point GMC, which lies in the excluded 
region. The region for the scenario of the gluino as the NLSP is specified by the hatched area, 
which is not accessible to this analysis. 
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spending to an integrated luminosity of 35 pb^ , collected with the CMS detector at the LHC, 
we find no evidence of such an excess. Our results are interpreted as 95% CL lower limits on 
GGM SUSY particle masses in a benchmark model. Additionally, the 95% CL upper limits on 
the cross sections of generic SUSY production as a function of squark / gluino and wino masses 
are also provided, allowing other SUSY scenarios to be constrained. These limits are the most 
stringent to date for a large range of parameters in this GGM SUSY scenario. 
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